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We study normal state properties of an interacting Fermi gas in an isotropic harmonic 
trap of arbitrary dimensions. We exactly calculate the first-order perturbation terms in 
the ground state energy and chemical potential, and obtain simple analytic expressions of 
the total energy and chemical potential. At zero temperature, we find that Thomas-Fermi 
approximation agrees well with exact results for any dimension even though system is dilute 
and small, i.e. when the Thomas-Fermi approximation is generally expected to fail. In the 
high temperature (classical) region, we find interaction energy decreases in proportion to 

T~2 J where T is temperature and d is dimension of the system. Effect of interaction in the 
ground state in two and three-dimensional systems is also discussed. 

§1. Introduction 

Recently Fermi degeneracy in neutral atomic gases has been attained for ^''K 
and ^Li^)"^). In these systems, one can control physical parameters with relative 
ease. Especially, the interaction strength can be varied in a wide range when Fes- 
hbach resonance is available^). It is thus expected that the exciting phenomena of 
Fermion paring either in momentum space (Cooper pairing) or real space (BEC of 
preformed Fermion pairs), leading to weak-coupling or strong-coupling superfluid- 
ity^)' This system is also characterized by the existence of trapping potential, 
which leads to discrete quantum eigenstates and shell structures. This makes ana- 
lytical calculations rather difficult. In many practical cases one uses Thomas- Fermi 
approximation for analytic calculations. This method is justified in the case when 
system has sufficiently large number of particles. Several properties ^'^^ ' such 
as single particle excitation and Cooper pair correlation is investigated in this ap- 
proximation. However it can't be always applicable when system is small and level 
discreteness play important roles. 

The purpose of this paper is to investigate the normal ground state properties of 
dilute interacting Fermi gas in a trap with use of the exact eigenstate of the harmonic 
trap without resorting to the Thomas-Fermi approximation. We give a systematic 
method to count up quantum number of harmonic potential when inter-particle 
interaction is contact type, and obtain a simple analytical expression of the total 
ground states energy and the chemical potential within the first-order perturbation 
theory. Our method gives exact energy shift caused by interaction in arbitrary 
number of particle as for the lowest order term in interaction strength. Moreover, 
this method makes it possible to investigate effect of interaction to structure of 



E-mail: hiroyuki@kh.phys.waseda.ac.jp 

typeset using VTpJ^.sty <ver.0.8> 



2 



H. Yoshimoto and S. Kurihara 



degenerate oscillator shell exactly. 



§2. GENERAL FORMALISM 



We consider a two-component Fermi gas with a short range interaction in an 
isotropic harmonic potential. The two components (hyperfine state) are labeled by 
a =11 and assumed to have equal concentration. We choose a unit system in which 
h = CO = m = 1, where m is the atomic mass, and uj is the frequency of an isotropic 
trapping potential. We start from the Hamiltonian H = Hq + .ffint; where 

H, = Y.j + ^r^)^a(r), (2-1) 

is the one-particle part, V'a(i') being the field operators of the component a. 

Hint =5/ V}(r)V|(r)Vi(r)^T(r), (2-2) 

is the two-body interaction Hamiltonian, and g is interaction strength. We don't use 
pseudopotential S(r)drr^^^'^^'> , and use a delta function for simplicity. Hq is easily 
diagonalized by writing the field operator as 

^a(r) = nE^7=i^«-, (2-3) 

i ni Y7r5n*!2"'' 

where d is dimension of the system, H^i{xi) is the Hermite polynomial, and n = 

(n^,,) is quantum number, for example n = {v} ,71^ ^n^) stand for d = 3, and 
[aani o^a'n'\ ~ ^oca'^nn'- Then Hq and i^int is written in the form 

d . 

Ho = J2J2i^'+2^^Laan, (2-4) 



Hint = 9^1 XI ^K'"2>'^3>'^4)«Tni«In2°i"3aTn4- (2-5) 



Here, ^(ni, n2, rza, 7x4) is a function of quantum numbers of trapped system and is 
given by 



u;(ni, 712, ns, 71.4) = (-1) 2 



'(2ni)!!(2n3)!! 



V (2n2)!!(2n4)!! 

foo f]U ,2 1.2 1.2 

I V^~"'^"'~"'^"''^~~^"'~"'^Y^^"3~"'(y)' (2-6) 

where L'!^{x) is the associated Laguerre polynomial. We can calculate this form using 
Fourier transformation Detailed properties of w{ni,n2, n^, n^) is investigated 
in Ref. 19). 
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§3. GROUND STATE PROPERTIES 

3.1. One- dimensional case 

We use the first-order perturbation theory, and obtain one-dimensional interac- 
tion energy for one-dimensional case in the form 

np 

Eintinp) = g ^ w{n,n,m,m). (3-1) 

n=0,m=0 

Here the function w{n,n,m,m) has complicated dependence on quantum numbers, 
showing considerable difference from momentum expansion ordinarily used to treat 
homogeneous system. We can nevertheless calculate Eq. (3-1) exactly as shown in 
Appendix A, with the result 

Above summation can be reduced further to an integral, yielding a simple form of 
energy 

EiN)^ — + ^m. (3-3) 

Here, N = 2N-^ = 2N\^ is total number of particles. Chemical potential is also 
calculated by the relation //(Aq,) = {E{Na + 1) - E{Na))/2, yielding 

1 2g 

- + A, + J 

2 vr^ 



f^{Na) = 77 + A, + -^Vm^. (3-4) 



Equations (3-3), (3-4) show that as A increases, ratio E\^^/Eq become small, thus 
one may see this system behaves like free particle system. However we should note 
the fact that even a weak interaction destroys Fermi-liquid nature in one dimension. 

3.2. Two-dimensional case 

In two-dimensional case, we have to deal with the degeneracy to calculate the 

energy. It turns out, however, as shown in Appendix B we can calculate the total 
energy exactly in this approximation when particle is filled up to Fermi level, with 
the result 

Einp) = \{nF + l){nF + 2){2nF + 3) 

+T|-(nF + l)(nF + 2)(2nF + 3). (3-5) 

The Fermi number up and total number A are related by np = ^{—3 -\- V4A + 1) 
therefore, energy is given by 

S(A) = ^V4A+T(1 + |-). (3-6) 
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Both the single-particle term and the interaction term have the same dependence on 
number of particle, similar to the results for homogeneous system. Mean value of 
chemical potential n{N) = is given by 

m(A^) = ^^^ ^i±^ (l + A). (3.7) 

In the following, we investigate the one particle energy shift. Let us suppose to add 
one more particle to system of which particles is filled up to Fermi level. Energy 
shift ^{hx, riy) is generally a function of degenerate level Ux-, Uy {hx + % = nj? + 1). 
But in this case, it doesn't depend on the coordinate, in other words, interaction 
between particles in the different energy levels doesn't lifts degeneracy, and is given 

by 

^(nj:,ny) = j^inp + 1). (3-8) 

Note that interacting term in Eq. (3-7) corresponds exactly to the above form. Next, 
we add one more particle having another hyperfine state to the system. Interaction 
between two particles on the same oscillator shell gives the deviation from Eq. (3-8). 
This energy shift S^{nx, Uy] rux, my) is given in Appendix C, with the result 

d^{nx,ny;mx,my) = 

-^^==(l + ilog4(ni. + l)) nx = mx = 
2^pf^(l + ilog4n,)(l + ilog4nj,) = / 
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T!-^{nx-mx) ' 



K(iJ^^^)KiiJ-^^^) Ux > rUx 



(3-9) 



Here mx,my (nix + niy = np + 1) are coordinates of the second particle, and K[x) 
is the complete elliptic integrals of the first kind. We find that the interaction be- 
tween two particles in the same oscillator shell lifts degeneracy mentioned above. 
IJ,{nx,ny;mx,my) have large value when two coordinate is close to each other, espe- 
cially it have maximum value when Ux = rUx = 01 Uy = my = 0. Eq. (3-8) and Eq. 
(3-9) imply two particles in the same shell make bound pairs in the same quantum 
number, when interaction is attractive and system is sufficiently dilute. 

3.3. Three-dimensional case 

Calculation can be done more or less similarly also in also three-dimensional 
case. The detail is given in Appendix D. The result for total energy is given by 

Einp) = \inF + l)(nF + 2f{nF + 3) 

g ^ / (2r + 3)!! \^ (2nF-2r + l)!! 
9(2^)1 (2r)!! ) {2n,-2r)}}. ' ^''^"^ 
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Equation. (3-11) can be also obtained from Thomas-Fermi approximation with 
weak interaction (| g |<C 1) i^)" ^3), 20) E^{N) = ^Jd^xn{x)^, where n{x) = 

— ^x^]^/^ is non-interacting particle density. It is possible to indicate the 
same procedure in one and two-dimensional case. As can be seen from Fig. 3, ex- 
act first-order perturbation calculation for the interaction energies are very close to 
the Thomas-Fermi approximation. Although the Thomas-Fermi approximation is 
ordinarily correct in the system, which has sufficiently large number of particles, 
our results show this approximation is applicable to any dimensional trapped sys- 
tem, even if the system is dilute as this approximation is generally expected to fail. 
Average chemical potential is given by 

128^2^ 



^^{N^) = (6iV„)i (1 + }^0{mS^y (3T2) 



We investigate energy shift more precisely as we have done on the two dimensional 
cases. One particle energy shift C{nx,ny,nz) is found to be 



V2g dJ'F 



47ri d^^F 



[p'^^''^FnMFnMFnM) 1/3=1 • (3-13) 



Where nx,ny,nz {rix + Uy + Uz = np + 1) is coordinate of degenerate states and 
Fn{l3) stands for Hyper geometric function F{—n, ^, 1; In this case, interaction 
lifts the degeneracy and Eq. (343) has maximum value when two of the quantum 
numbers are zero. The maximum value is 

CK, Uy, n,) - ^ (i-nf + ^ni) . (3-14) 

This form reproduces results of Eq. (3T3) have a minimum value at = 

Hy = Uz, namely one particle wave function assumes spherically symmetric shape. 
We calculate the minimum value numerically and find it to be roughly 0.89 times 
lower than maximum value. 

§4. FINITE TEMPERATURE 

At the finite temperature, number of particles and energy of d-dimensional sys- 
tems is written in the form 

oc 

iV = 2 ^ D(n)/(n) (4-1) 

n=0 

00 J 

E = 2j2D(.n){n+-)f{n) 

n=0 

d 00 

+^11 E Hni,mi)f{n)f{m), (4-2) 

i=l ni=0 
mi=0 
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here, /(n) is distribution function, n = ni + ..n^, and D{n) = J2ni+..nd=n ^ density 
of state. We restrict our discussion to high temperature, i.e. well above the Fermi 
degeneracy. Then, we use Boltzmann distribution function and calculate Eq. (4-1) 
and find to be 



N = 2 



(1-e t] 



^ 2r'^e^('*-i). (4-3) 

Here, we put on Boltzmann constant /eg = 1 for simplicity. Using the Eq. (4-3), we 
can calculate Eq. (4-2) which is given in Appendix E with the result 



79T + — tanh2 (- 

2tanh(|) 4(27r)i 



^dNT+— ^. (4-4) 

4(47rr)i 

We have thus obtained both for T = and T ^ oo limits. These results may be 
useful in constructing an interpolation formula for finite temperatures. 



§5. SUMMARY 



In summary, we have studied an interacting Fermi gas in arbitrary dimensional 
isotropic harmonic trap. We have used first-order perturbation theory. In contrast to 
the homogeneous system, first-order perturbation term has non-trivial contribution 
to the system. Dependence on the number of particles, dimension, and temperature 
is characteristic of the trapped system. We have calculated energy and chemical 
potential and found that Thomas-Fermi approximation is applicable when the sys- 
tem is close to the ground state, even if the system is dilute and small, i.e. when 
the Thomas-Fermi approximation is generally expected to fail, and we have found 
interaction energy decreases in proportion to T~'^ , when the system is well above 
the Fermi degeneracy. 
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Appendix A 

Derivation of Eq. (3-3), Eq. (3-4) 



One easily find Eq. (3-1) is written in the form 

np np 



Ei^tinp) = / dkk-^e-'' V Lnik) V L^(fe) (A-1 

V27r Jo ^0 ^0 
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We use the relation 



^^(-) = E^S7^i?^t.(-), (A-2) 



- riria-b) 

here, r{r) is gamma function, we put on a = 1 and 6 = and find to be 

Ei^iinp) = r dkk-\e-^Ll^{kf. (A-3) 
V27r Jo 

We use Eq. (A-2) again, put on a = 1 and h= \ and use the relation 

r dxe--x-Ll{x)Ll,{x) = + ^ + (A-4) 
Jo 

One easily finds Eq. (3-2) from Eq.(A-4). We use Wallis's formula (2r- l)!!/(2r)!! ^ 
Xj^fwr and replace the summation with integral, then interaction term of Eq. (3-2) 
is reduced to be 

^ ' dr ^ (A-5) 



\/27r2 Jo ^np - r' 

We obtain Eq. (3-3) directly from (A5). Interaction term of chemical potential is 
written in the form 

/^int(iVa) = (^int(iVa + 1) " £^int(iVa))/2 

= ff X] f{n,ri,nF,nF). (A-6) 

n=0 

We use the relation Eq. (A-2) and Eq. (A-4) and find to be 

g "^V(2r-1)!!\2 



9 f{2r-l)U\ 
(2nF + l-2r)!! 



{2nF + 2 - 2r)! 



(2r + 1) 



' dr ^ (A-7) 



vr^ 7o + 1 - r 

We obtain Eq. (3-4) directly from (A-7). 

Appendix B 

Derivation of two-dimensional energy 



Two-dimensional interaction energy is written in the form 



Eint{N)=g^J2 E h{nx,m^)h{ny,my). (B-1) 



n=0 m=0 Ux +ny =n 
mx+my=m 
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Here, h{n,m) represents w{n,n,m,m.). Wc use the relation 

^ Lr^,{xi)..L„^{xd) = Lt\xi + .. + Xd). (B-2) 

ni+..+nd=n 

We put on d = 2, then Eq. (B-1) is reduced to be 

np np 



„ n ™ n J — OO 



n=0 m=0 ■ 



We replace above Cartesian coordinates valuable with polar coordinates valu- 

ables (r, ^) and use Eq. (A-2), then Eq. (B-3) is reduced to be 

We apply Eq. (A-2) and Eq. (A-4) to above form, and we obtain Eq. (3-5). Note 
that three-dimensional energy can be also calculated as well. 

Appendix C 

Derivation of Eq. (3-9) 



Eq. (3-9) is written in the form 

SC{nx,ny;mx,my) = 

gh{nx, my)h{nF -M - n^, -M - my). (C-1) 

We use Eq. (A-2) and Eq. (A-4) and find 

" (2r-l)!! 
/i(n,m) = V27r_^ (2r)!! 

(2m - 2n + 2r - 1)!! (2n - 2r - 1)!! 
{2m-2n + 2r)\\ (2n-2r)!! 

If n = m above form is approximated to 



(C-2) 



/I 1 
h(n,n) ^ \/2^ ^= + ^ / dr 



1 r^/n — r 
^^ + ilog(4n)Y (C-3) 



n V TT 



else if m > n approximated to 

h{n, m) ^ V2^— / tir 

7r2 Jo vr(n — rj(m — n + rj 

2^2 



TTx/m — n \ m — n 
We use Eq. (C-3), (C-4) and obtain Eq. (3-9). 



(C-4) 
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Appendix D 

Derivation of Eq. (3-13), (3-14) 

We consider adding one more particle to three-dimensional system of which 
particles is filed up to Fermi level. Then, energy shift is 



^ ^ h{mx,nx)h{my,ny)h{r 



^ ^ m=0 mx+my-\-mz=m 

OO n—— 2 

= ^ / drrV^'^' / ' # / ' d^sin^L^ (^) 

TT'^ Jo Jo Jo 2 

2 2 

X -^nx (y sin^ ^ cos^ v)Lny (y sin^ 6* sin^ ip) 

2 

xL„,(^cos2^). (D-1) 
We expand the Laguerre polynomials and calculate angular integrals 

2 



/ dip i (19 sin 6 Ln^(— sin^ 9 cos^ (p) 
Jo Jo 2 

xL„^(— sin^ ^sin^ cos^ 9) 

, "y riz / X / x / 

-1 / 1 \i / \ / \ / "-^ 



^ ii=0Z2=0«3=0 ^ X / \ y / \ z 



here, ^ ^ ^ is binomial coefficient, and we puts on I = + ly + Iz- Radial integral 



I 

is also calculated and find to be 



r drrl^^e-^Ll^{r) = ^^^,"^ + "f r(^ + Z). (0-3) 
Jo ni7!r(| — /) 2 



Prom Eq. (D-2) and Eq. (D-3), we find 



47r up. 



(D.4) 
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here, b{n,l) is (-1)^( ^ ) r(/ + i)//!. Eq. (3-13) is derived from above form. When 
Uy = Hz = 0, Eq. (D-4) is reduced to be 

+ 1, 0, 0) = {n, + 1)^ "g' fci^ 

(2/-5)!!(2nf-2Z + 3)!! 
(20!! (2nir + 2-2Z)!! ^ ^ 

We find that ^^^^T)^]'' *'^(27)l!'' decreases in proportion to on the other hand (2nF+2^-^0!! 
decreases in propotion to {np—l)^ and have nearly constant value if I « np- Hence, 
we find 

X . ^V2g "p.^ {2l-iy.l 

C(n^ + 1,0,0) -(n^ + l)|^ g 

x^^|^(n| + 0(n7)) 
4a/2 3 1 

= 9^»F + 0(n^). (D-6) 



Here, wc use the relation TrJ^f^o ^"^[21)11' ^'^(2i)n ' ~ ^/^' ^® obtain the first term of 
Eq. (314). Coefficient of second term is approximated by summation of Eq.(D-5) 
up to i = 2. 

Appendix E 

Derivation of Eq. (4'4) 



At the high temperature i.e. Boltzmann distribution is applicable, d-dimensional 
energy is given by 

00 7 
Eo = 2j2Din){n+-) (E-1) 

n=0 

(l_e-T)<i+i 
^de^{f^-i)(^l + e-^^T'^+\ (E.2) 

Kinetic term of Eq. (4-4) is obtained from Eq. (E-2) and Eq. (4-3). On the other 
hand interaction term is given by 

00 00 

^int = 7^ E E E e-^(-+--'^+i) 

V / n=Om=0 n\+..nd=n 
mi + ..m^=m 



xe-^("*i+--"^d-M+f) P dxi..dx<ie-5(^?+-^d) 
J —00 
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x'^ x^ x^ x^ 

^-^mi (-2-)--^'md (-2-)-^ni (-2-)--^'nd (-2")- 



(E-3) 



We use the relation (B-2) then, above form is reduced to be 

g 



int 



(27r)'^ r(f + 1) 



oo 2 



drr" 



.2 \ 2 



n=0 

We use Eq. (4-3) and the relation 

oo _ 



(E-4) 



xt 

e 1-* 



and find to be 



int 



(27r)'^ r(f + 
We obtain Eq. (4-4) from Eq. (E-6). 



2 

fOO L 

— A^^ / drr'^~'^e ^tanM^j^) 
1) Jo 



(E-5) 



(E-6) 
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Fig. 1. Energy shift caused by interaction between two particles in the same energy level np + l = 41 
for g = 1. {nx,ny) = {nx,nF + 1 — rix) and {mx,my) = {mx,nF + 1 — rrix) stand for two- 
dimensional degenerate states of the two particles. 




Fig. 2. Energy shift caused by interaction between two particles on the same two-dimensional 
degenerate states (nx,ny) = {nx,nF + 1 — rix) for g=l. 
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Fig. 3. Ratio as a function of Fermi level np for the interaction energy using Thomas- Fermi ap- 
proximation E^{nF) and exact first-order perturbation, d represents dimension of the system. 




Fig. 4. Single particle energy shift of the three-dimensional system caused by inter-shell interaction 
for g — 1. {nx,ny,nz) = {nx,ny,nF + 1 — rix — Uy) is degenerate state of the particle. 



